Resummation for W and Z production at large px 
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Soft-Collinear Effective theory is used to perform threshold resummation for W and Z production 
at large transverse momentum to next-to-next-to-leading logarithmic accuracy including matching 
to next-to-leading fixed-order results. The results agree very well with data from the Tevatron, 
and predictions are made for the high-pT spectra at the LHC. While the higher-log terms are of 
moderate size, their inclusion leads to a substantial reduction of the perturbative uncertainty. With 
these improvements, the PDF uncertainties now dominate the error on the predicted cross section. 



The production of an electroweak boson at large 
transverse momentum is arguably the most basic hard- 
scattering process at hadron colliders. In fact, it was one 
of the first for which the next-to-leading order (NLO) 
perturbative corrections were computed [l|-|3|- By now 
the complete a 2 corrections to vector boson production 
are known, but since the pr-spectrum starts at 0(a s ) 
its theoretical accuracy has not been improved since the 
nineteen eighties. Given that Z's and W's at high trans- 
verse momentum provide an important background to 
new physics searches and a standard way to calibrate jet 
energy scales, it is important to have good theoretical 
control of their cross sections. 

In the absence of a full NNLO computation, a way to 
improve predictions is to compute the contributions to 
the cross section which arise near the partonic threshold 
from the emission of soft and collincar gluons. For very 
large p^, these corrections become dominant and must be 
resummed to all orders in perturbation theory. However, 
even away from this region, the corresponding terms yield 
in many cases a good approximation to the full cross 
section. For vector boson production, this resummation 
has been performed to next-to-leading logarithmic (NLL) 
accuracy m 

In the present paper, we use Soft-Collinear Effective 
Theory (SCET) [7H9[ to perform the resummation of the 
threshold terms to NNLL accuracy. A detailed deriva- 
tion of the factorization theorem and phcnomenological 
analysis of the closely related process of direct photon 
production were presented by the authors in [l(J. This 
paper is a generalization of those results to include vector 
boson masses which significantly complicates the calcu- 
lations. 

In the partonic threshold region, where the final state 
has an invariant mass much lower than the transverse 
momentum, the cross section for a given partonic initial 
state / to produce an electroweak boson V factorizes as 



do - / 
dsdi 



af(s,i)H I (s,i,M v , f i) 



where the partonic Mandelstam variables are s = (p± + 
P2) and t = (p\ — q) 2 , with q the vector boson momen- 
tum, with q 2 = My. We have factored out the Born 
level cross section af(s, t). The parton momenta p^ and 
P2 carry momentum fractions X\ and X2 and are related 
to the hadron momenta via p^ = x\P^ and p^ — X2P2 ''■ 
The hadronic cross section is obtained after convoluting 
with parton distribution functions (PDFs) and summing 
over all partonic channels. 

At the partonic threshold the mass of the final-state 
jet m 2 x = (pi + P2 — q) 2 vanishes. Since the full final 
state also involves the remnants of the scattered hadrons, 
this condition does not imply that the hadronic invariant 
mass M| = (P1+P2— q) 2 vanishes, unless the momentum 
fractions x\ and X2 are close to unity. Because there is 
no phase space for additional hard emissions near the 
threshold, the only corrections to the cross section arise 
from virtual corrections, encoded in the hard function 
Hr(s,t, My, A*), and soft and collinear emissions given 
by Si(k,fi) and Ji(p 2 ). The convolution over the soft 
momentum in (fT|) arises because the partonic jet mass is 
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x dkJ I (m 2 x -2E J k)S I (k,Li), (1) 



where Pj and kg are the collinear and soft momenta in 
the jet, Ej is the jet energy and k = pj-ks/Ej. Since the 
reaction in the threshold region proceeds via Born-level 
kinematics, only two partonic channels are relevant: the 
Compton channel qg — > qV and the annihilation channel 
qq -> gV. 

A detailed derivation of the factorization formula (fTJ) 
was given in [loj . That paper focused on photon pro- 
duction, but exactly the same jet and soft functions are 
relevant also for W and Z production. Explicit expres- 
sions for these functions in both the Compton and an- 
nihilation channel can be found in (Io| . together with 
the relevant anomalous dimensions. The non-zero boson 
mass only enters the hard function Hj(s,t, My, fj) and 
modifies the kinematics. The hard function is given by 
the virtual corrections to the corresponding hard scatter- 
ing channel. To obtain the function at NLO, one needs 
the interference of the one-loop amplitude with the tree- 
level result, which is given in equations (A. 7) to (A. 9) of 
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FIG. 1: Resummation by RG evolution. 

0. These expressions correspond to the bare result for 
the hard function. After rcnormalization, we obtain 



Hqq (s, i, M v , n) = 1 + -p 



-in 



s 



21n^ ( CUln£ + 3Cf] -6C F \n-^ 
s \ tu J Mp 
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A£(M)/4 (3) 



for the hard function in the annihilation channel, with 



, t K 



The extra piece AC(s,t) is the part of the virtual cor- 
rections not proportional to the Born cross section. It 
is finite and given in the last five lines of (A. 9) of Q. 1 
The virtual corrections for the Compton channel are re- 
lated to the above by crossing, see [2| for the necessary 
relations. 

Instead of integrating over the momentum fractions X\ 
and x 2 of the two partons involved in the hard scattering, 
we follow [l| and change variables to 



<Lx\&x 2 6(rn x ) 



dxi 



jJn X\S + U - M$r J 



dm 2 x , 



with U = (P 2 - q) 2 = M 2 - yfs^/M 2 + p 2 ^ and 

m max = U + X 1 (M| - U) , X mh 



(5) 



M\-u 



In the above variables, the expansion around partonic 
threshold m x — is performed at fixed Xi, which is 
the same as fixed t = (x\P\ — q) 2 . However, simply 



1 The result is also given in Q| but the sign of the AC terms in 
expression (A4) in this reference is incorrect. We thank R. Gon- 
salves for helping us to resolve this discrepancy. 
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FIG. 2: Scale setting. The plots are for W + bosons, but the 
qualitative features are the same for all bosons. 



using ([5]) would be problematic, since the choice of ex- 
pansion variables is not invariant under crossing p\ <-> p 2 , 
which would induce unphysical asymmetries in the rapid- 
ity spectrum. To avoid this, we symmetrize by perform- 
ing the expansion around threshold twice: once at fixed 
x\ using ([5]) and once at fixed x 2 ■ Our result for the cross 
section is the average of the two expansions. 

To perform the resummation of the soft and collinear 
emissions, the hard, jet and soft functions must be eval- 
uated at their own characteristic scales and then evolved 
to the factorization scale using the renormalization group 
(RG) where they are combined with the PDFs. This is 
illustrated in Figure [1] The running is straightforward 
using the Laplace space formalism developed in 

The most naive approach to scale setting would be to 
set the jet scale fij equal to the partonic jet mass mx- 
However the partonic jet mass is not an observable; it 
is integrated over in the convolution with the PDFs. So 
setting fij = mx one will encounter the Landau pole 
in the strong coupling constant, since the partonic jet 
mass mx can become arbitrarily small. The problematic 
choice fij = mx is inherent in the traditional formalism 
for resummation. To avoid the Landau pole previous 
work did not perform the resummation to all orders and 
instead just computed the singular terms in the partonic 
cross section at the next, or the next two orders in per- 
turbation theory 0, 0. In our work, we resum to all 
orders but choose [ij to be the average-jet mass. Since 
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the partonic cross section is convoluted with the PDFs, 
the average jet mass must be calculated numerically. 

To obtain the proper scale for each ingredient of the 
factorization formula Q]), we use the numerical procedure 
advocated in [l2| • We evaluate the factorization theorem 
(fTJ) numerically at a fixed scale [i and study individually 
the impact of the NLO corrections from the hard, jet 
and soft functions. The scale variations are shown on 
the top panel in Figure [5] Note that the jet and hard 
function variations have natural extrema. These extrema 
are shown as the points in the lower panel for the jet scale. 
The solid curves are a reasonable approximation to these 
points, given by 
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which we use instead of the exact extrema for simplicity. 
We also set /i s = Hj/fJ,h, as dictated by the factorization 
theorem. By choosing scales close to these extrema, we 
minimize the scale uncertainty. 

The scale setting procedure beautifully illustrates the 
power of the effective field theory approach. In a fixed- 
order computation, the hard, jet and soft corrections can- 
not be separated and are included at a common value 
of the renormalization scale. This is shown in the NLO 
curves in the top of Figurc[2J In this case the scale depen- 
dence is monotonous. Because there are multiple relevant 
scales in the problem, there is no natural scale choice at 
fixed order. But there are natural choices when (i is split 
into hard, jet and soft. For illustration, we also show 
in the bottom panel of Figure [5] the popular scale choice 
fij = \/p'x + My, which is a bad fit to the jet scale. 

The most natural choice for the factorization scale (if 
would be at or below fi s , since this scale defines the 
boundary between the perturbative and non-perturbativc 
part of the process. However, since all PDF fits were per- 
formed with fi f set equal to the hardest scale in the pro- 
cess, we will follow this convention and use /if = /x/, as 
our default value. Note that this implies that we use the 
RG to run the jet and soft functions from lower to higher 
values, in contrast to the situation depicted in Figure [TJ 

In order for our results to contain the full NLO cross 
section we match the resummed result to fixed order. 
The matching is straightforward since the resummation 
switches itself off when we set all scales equal, (ih = fij = 
fi s = (if. Doing so in the NNLL result yields all loga- 
rithmically enhanced terms at one loop level, which we 
denote by NLO S i ng . The matched result is obtained by 
adding the difference between the full NLO result and the 
singular terms NLO s ing to the NNLL resummed result. 
We denote the matched result by NNLL+NLO. To com- 
pute the NLO fixed-order result, we use the code QT [l3j ]. 
and have verified that it agrees with mcfm [l4[ . 
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FIG. 3: Comparison to DZero results [Is| |. 
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FIG. 4: Prediction for the combined W + and W cross sec- 
tions the LHC (7 TeV). 



Most of the ingredients necessary to go to N 3 LL ac- 
curacy are already known. The 3-loop anomalous di- 
mensions are all known. We use these, along with the 
two- loop jet function constants [H, [lj| and the Pade ap- 
proximant for the 4-loop cusp anomalous dimension, to 
get our most accurate prediction, which we denote by 

N 3 LL par tial. 

For fixed order, we set the renormalization and factor- 
ization scales equal to fj,h- Uncertainties are estimated 
by varying by a factor of 2 around the defaults values 
([6]) and extracting the maximum and minimum values. 
For the final scale variation error bands on the resummed 
distributions, we add the jet, hard, soft and factorization 
scale uncertainties in quadrature. All numerical predic- 
tions arc computed using MSTW2008NNLO PDFs [l7j|, 
with a s (Mz) = 0.1171. For the electroweak parameters, 
we use a = 1/127.92, sin 6 W = 0.2263, M w = 80.40 GeV, 
M z = 91.19 GeV. 

In Figure [3l we show the pr spectrum of the Z-boson 
at the Tevatron in comparison to results of the DO exper- 
iment [Til ] . Our results agree well with the experimental 
results but have significantly smaller uncertainties, in the 
region of high transverse momentum. 



4 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



LHC (7TEV) 




N 3 LL p „ llal +NLO 






<x (pp ^ Z) 
oMpp -> W*) 


PDF uncertaiiitie^^^rrr^^^^ 






o-(pp -> W-) 




o-(pp ^W+) 



Tevatron 
W ± Z 



LHC at 7TeV 
W ± Z 



200 400 600 800 1000 1200 

Pt 



FIG. 5: Prediction for ratios of cross sections for different 
gauge bosons at the LHC (7 TeV). 



In Figure |H we give results for the production of W 
bosons at the LHC. The perturbative uncertainty on our 
result is much smaller than the PDF uncertainty, which 
implies that these results could be used to obtain more 
precise determinations of the PDFs, once the experimen- 
tal results become available. 

In Figure [5j we give the ratio of the W~ to W + cross 
section, as well as the ratio of Z's to W's. Most theoret- 
ical uncertainties drop out in these ratios, making them 
precision observables. Roughly half as many W~ as W + 
bosons are produced at the LHC, which is just a reflection 
of the quark content of the protons, while the Z to 
ratio is essentially given by the ratio of the electroweak 
couplings. Note that the PDF uncertainty is larger on 
the ratio involving the W~ , since the W~ is sensitive to 
sea-quark PDFs. 

In Table U we give numerical results for integrated 
cross section with pt > 200 GcV at the Tevatron and 
LHC. In this table, we have included rows NLO S i ng . and 
NNLOsing. which refer to the resummed distributions ex- 
panded to fixed order. Indeed, Figure [2] and the pre- 
scription © make it clear that there is not a very large 
hierarchy between the different scales. For example, at 
the Tevatron, for pt = 200 GeV, the ratio between the 
scales is fih/fJ-s ~ 5. Since the logarithms of the differ- 
ent scales are of moderate size, one can re-expand the 
resummed result in a s and will find that most of the 
effect of the resummation is captured by the singular 
terms. NLO s ing. refers to the expansion of the NNLL 
result (without matching) to a 2 s . NNLO S m g . refers to 
the expansion of the N 3 LL par ti a i result to o? s . Compar- 
ing NLO to NLO S i n g. we see that for Z-production with 
p T > 200GeV at the Tevatron (LHC), we find that 82% 
(70%) of the NLO correction to the cross section is due 
to soft and collinear gluon emission. Thus we expect that 
SCET reproduces most of the NNLO result at high px- 

Finally, there is one important effect that we have ne- 
glected, namely electroweak Sudakov logarithms. These 
can be large and must be included. Already at the Teva- 
tron these are non-negligible and lower the cross section 
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TABLE I: Cross section o(j>t > 200 GeV) (in picobarn) using 
different approximations, see text. In addition to the scale un- 
certainties shown in the table, there is a relative PDF uncer- 
tainty of 3% (2%) for W ± production at the Tevatron (LHC) 
and 5% (2%) for Z production . 



at p T = 300 GcV by about 10%. At the LHC at 14TeV, 
the effect is about —15% at pt = 500 GcV, and around 
-30% for pt = 2 TeV To get a really accurate 

comparison to data, the electroweak corrections must be 
included. 

It would be interesting to obtain NNLL predictions for 
less inclusive quantities. In particular, to compare to the 
preliminary CMS results for the Z boson pt spectrum 
[19(, based on 36 pb -1 , one would need results which are 
differential in the lepton momenta to account for the ex- 
perimental cuts. To obtain these results, one has to com- 
pute the hard function for arbitrary boson polarizations, 
which is in progress [22| . Putting in cuts on the hadronic 
side is more complicated, but would be quite interest- 
ing, since it would allow to compute the production in 
association with a jet. 

The LHC has just delivered its first inverse femtobarn 
of data at 7 TeV. Comparing our results to this data 
will allow for precision tests of the standard model at the 
highest energies ever produced in a collider. 
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